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Abstract. The theoretical prediction of micro- and macroturbulence and <fn,ac) as a 
function of stellar parameters can be useful for spectroscopic work based on ID model at- 
mospheres in cases where an empirical determination of ^^^^ is impossible due to a lack 
of suitable lines and/or macroturbulence and rotational line broadening are difficult to sepa- 
rate. In an effort to exploit the CIFIST 3D model atmosphere grid for deriving the theoretical 
dependence of fn,ic and ^mac on effective temperature, gravity, and metallicity, we discuss 
different methods to derive ^^ic from the numerical simulations, and report first results for 
the Sun and Procyon. In both cases the preliminary analysis indicates that the microturbu- 
lence found in the simulations is significantly lower than in the real stellar atmospheres. 
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1. Introduction 

3D hydrodynamical simulations of stellar 
surface convection provide a physically 
self-consistent description of the non-thermal 
velocity field generated by convection, 
overshoot, and waves. This is one of the 
great advantages over classical ID model 
atmospheres where the properties of the 
photospheric velocity field need to be specified 
empirically in terms of the free parameters 
^mic and ^mac- Evcu if these parameters are 
irrelevant in the context of 3D hydrodynam- 
ical model atmospheres, the latter may be 
used to predict the magnitude of and 
^mac as a function of the stellar parameters 
(^eff, logg, [M/H]). The CIFIST 3D model 
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atmosphere grid (Ludwig et al., this volume) 
provides a suitable database for this purpose. 
In the following, we investigate different 
methods to derive the parameter ^mic (and 
subsequently ^mac) from our 3D model at- 
mospheres which were comp uted with the 
CO^BOLD codJl (jj evtag, St effen. & DorchI 
l2002t IWedemever et al.i i200|). The Sun 
(reff=5780 K, logg=4.44, [M/H]=0) and 
Procyon (reff=6500 K, logg=4.0, [M/H]=0) 
serve as benchmarks for the present study. 

2. Methods to derive ^mic and ^mac 

All methods described below presume that the 
photospheric velocity field may be character- 
ized by a single, depth-independent value of 

' http://www.astro.uu.se/~bf/co5bold_main.html 
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^mic and ^mac, and that both micro- and macro- 
turbulence have an isotropic Gaussian proba- 
bility distribution of the line-of-sight velocity, 
P{v) ~ e\p(-v^/^^). Synthetic spectral hues 
serve as a diagnostic tool to probe ^mic and ^mac 
via the total absorption (equivalent width W) 
and the shape of the line, respectively. 

Method 1 (Ml) is considered the most 
accurate procedure to extract the microturbu- 
lence parameter from a 3D numerical convec- 
tion simulation. Unlike the other methods de- 
scribed below, it relies only on the 3D model, 
and yields a value of ^mic for any individual 
spectral line, thus allowing to map ^midW). 
Given the spectral line parameters, the line pro- 
file is computed from the 3D model with dif- 
ferent velocity fields: (i) using the original 3D 
hydrodynamical velocity field, and (ii) replac- 
ing the 3D hydrodynamical velocity field by an 
isotropic, depth-independent microturbulence, 
like in classical ID spectrum synthesis, but re- 
taining the full 3D thermodynamic structure. 
Now the microturbulence associated with the 
considered spectral line, is defined by 
the requirement of matching line strengths: 
W3oihydro) = W'sdC^^,), where l^aDCliydro) 
and WsDC^mic) are the equivalent widths ob- 
tained in steps (i) and (ii), respectively. 

Once ^mic is determined as explained 
above, the macroturbulence associated with the 
considered spectral line is defined by minimiz- 
ing the mean square difference between the 3D 
line profile computed in step (i) with the hydro- 
dynamical velocity field, and the modified 3D 
line profile obtained in step (ii) with ^J^j^ and 
subsequent convolution with the macroturbu- 
lence velocity dispersion (fmac)- The IDL pro- 
cedure MPFTT is used to find the solution 
that gives the best fit to the original 3D profile. 

Method 2a/b (M2a/h) The idea of this 
method is to replace the modified 3D models 
(furnished with a classical microturbulence ve- 
locity field) used in step (ii) of Method 1 with 
an arbitrary ID model atmosphere. However, 
this concept does not work for a single spec- 
tral line, but instead has to rely on a set of 
spectral lines ranging from weak to partly sat- 
urated. In Method 2 the set of Unes is gen- 
erated from a curve-of- growth, i.e. all lines 



share the same atomic parameters except for 
the /-value, which controls the line strength. 
Given the set of fictitious spectral lines and 
a ID model atmosphere, we first compute 
for each line i the 1D-3D abundance difl'er- 
ence A log €i{ f mic), defined by the condition 
Wzui^oggfi) = WiD(^mic,logSf/;- + Aloge,). 
i.e. Aloge, is the logarithmic abundance dif- 
ference between the abundance derived from 
the ID model by fitting the equivalent of the 
3D profile, and the true abundance assumed 
in computing the 3D spectral line. In M2a 
we compute the slope of the linear regression 
to the data set [W^u, Alog ei(^niic)K and de- 
fine by the condition of vanishing slope, 
I3{ ^ Q Alternatively, in M2b the micro- 
turbulence is required to produce the minimum 
standard deviation of {A log e,}, a{ ^*') = min. 
Once ^mic is determined, the macroturbulence 
associated with each individual spectral line, 
^ac' is found by fitting the original 3D line 
profile with the corresponding ID line profile 
(with fixed ^V'^) through variation of fmao m 
essentially the same way as in Method 1. 

Method 3a/b (M3a/b) is very similar to 
Method 2, except for utilizing a sample of real 
spectral lines of different strength (and differ- 
ent wavelength, excitation potential, etc.) in- 
stead of a set of fictitious lines lying on a single 
curve-of-growth. Adjusting the value of ^mic 
to minimize the difference in A log e, between 
weak and strong lines obtained from the pre- 
ferred ID model results in (zero slope 
condition). Similarly, minimizing the overall 
dispersion of Aloge,- gives ^^l- Method 3 
corresponds to the classical definition of ^mic- 
Again ^[^^^c is found as in Method 2. 

Note that all 3 methods have the advantage 
that errors in log gf cancel out. Only Method 3 
can also be applied to observed stellar spectra; 
in this case the accuracy of the log gf values 
is crucial. Obviously, the results depend on the 
selected spectral lines, and, for Methods 2 and 
3, on the choice of the ID model atmosphere. 

In principle, it should be possible to derive 
^mic and ^mac directly from evaluating the 3D 
hydrodynamical velocity field without resort- 
ing to synthetic spectral Unes. But a suitable 
procedure (Method 4) has yet to be developed. 
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3. First results: Sun and Procyon 



Figures [T] and |2] show the determination of the 
microturbulence parameter from a 3D model 
atmosphere of the Sun and Procyon, respec- 
tively, according to the three methods de- 
scribed above. We considered 5 sets of ficti- 
tious iron lines (Fe 1 0, and 5 eV, Fe II 0, 5, and 
10 eV, used with Ml, M2), and one sample of 
real Fe I and Fe II lines each (used with Ml, 
M3). The most obvious result is that the mi- 
croturbulence derived from the flux spectra is 
systematically higher than that obtained from 
the intensity spectra, in agreement with obser - 
vational evidence (e.g. iHolweger et al.l[T978l) . 
Moreover, the derived value of ^mic depends 
on the type and strength of the considered 
spectral line. Method 1 clearly reveals that 
high-excitation lines 'feel' a lower microtur- 
bulence than low-excitation lines. The depen- 
dence of ^inic on line strength is non-trivial (see 
Fig.ffl lower panel). In general. Methods 2a 
and 2b give very similar results, which more- 
over agree well with the results obtained from 
Method 1 for the stronger lines on the curve- 
of-growth. In contrast. Methods 3a and 3b can 
give rather discordant answers, if the sample of 
spectral lines is not sufficiently homogeneous 
(see again lower panel of Fig.[T]l. 



In Table[T] the theoretical predictions of 
micro- and macroturbulence based on our hy- 
drodynamical model atmospheres for the Sun 
and Procyon are summarized and confronted 
with empirical results from the literature. This 
preliminary analysis indicates that the theoret- 
ical predictions of ^mic fall significantly below 
the classical empirical estimates, for both so- 
lar intensity and flux spectra, and even more 
clearly for Procyon, as demonstrated unmistak- 
ably in Fig.|3] On the other hand, the macrotur- 
bulence derived from the CO^BOLD models is 
somewhat larger than deduced from observa- 
tions, such that the total non-thermal rms ve- 
locity fturb (columns (6) and (7) of Tab.[T]i, and 
hence the total line broadening, is very similar 
in simulations and observations. 




7.4 ^ 



2 4 6 8 10 12 

WJpm] 

Fig. 3. Iron abundance derived from individual 
Fe II lines using a 3D CO' BOLD model with 
reff=6500 K, log g=4.0, [M/H]=0, with log gf 
values taken from Melendez & Barbuy (2009), and 
va n der Waals l ine broadening treated according 
to lBarklem et al.l ( [T998) . The strong trend of abun- 
dance with equivalent width indicates an apparent 
lack of small-scale turbulence in the hydrodynami- 
cal velocity field of the Procyon model. 



4. Conclusions 

We have developed different methods to extract 
the parameters and ^mac from 3D hydrody- 
namical simulations, and found all methods to 
give consistent results. As expected, ^mic and 
^mac depend systematically on the properties 
of the selected spectral lines. Our preliminary 
analysis based on state-of-the-art CO^BOLD 
models for the Sun and Procyon indicates that 
the microturbulence seen in the simulations 
is significantly lower than measured in the 
real stars. This suggests that the velocity field 
provided by the current 3D hydrodynamical 
models is less 'turbulent' than it is in real- 
ity. It seems that the simulations predict essen- 
tially the correct total rms velocity, while un- 
derestimating the small-scale and overestimat- 
ing the large-scale fluctuations. If confirmed, 
this deficiency implies a systematic overesti- 
mation of 3D abundances from stronger lines. 
Apparently, our findings are in conflict with the 
conclusions by lAsplund et al.l (l2000l) . Further 
investigations are necessary to find the basic 
cause of this disagreement and to clarify possi- 
ble implications for high-precision abundance 
studies. 
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gt57g44n58_35x35: disk-integrated spectrum 
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Fig. 1. Overview of ^n,i;. derived from our standard 3D solar model atmosphere (gt57g44n58) determined 
with Ml, M2, and M3 for a total of 60 different Fe I (filled circles) and Fe II (plus signs) lines. Connected 
symbols represent fictitious lines sharing the same curve-of-growth, individual symbols indicate real spec- 
tral lines. The results of M2a,b (red, blue) and M3a,b (green) are shown to the right of the vertical line 
in each panel for two different ID models (smaller symbols: (3D), larger symbols: HM). Top and bottom 
panels refer to disk-center and full-disk spectra, respectively. 
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d3t65g40mm00n01_35x35: disk-center spectrum 
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d3t65g40mm00n01_35x35: disk-integrated spectrum 
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Fig. 2. As Fig.[Tl but for a 3D model atmosphere representative of Procyon, using a total of 56 different Fe I 
(filled circles) and Fe II (plus signs) lines. ID models used with M2, M3: (3D) (smaller symbols) and LHD 
[Q'mlt=0.5] (larger symbols). Top and bottom panels refer to disk-center and full-disk spectra, respectively. 
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Table 1. Empirical values of ^^ic and ^mac from the literature, compared with the theoretical 
results derived in this work from 3D CO^BOLD hydrodynamical model atmospheres for the Sun 
and Procyon (weighted average over Methods 1, 2 and 3). 



Atmosphere / 




[km/s] 




[km/s] 




f^Lc)/2 [km/s] 


Model 


disk-center 


full-disk 


disk-center 


full-disk 


disk-center 


full-disk 


Sun, observed" 


1.00 ±0.15 


1.35 ±0.15 


1.63 ±0.15 


1.90 ±0.15 


1.35 ±0.10 


1.65 ±0.10 


3D solar model 


0.70 ±0.10 


0.95 ±0.15 


1.85 ±0.45 


2.30 ± 0.30 


1.40 ±0.30 


1.75 ±0.20 


(gt57g44n58) 














Procyon, observed* 




2.10 ±0.30 




4.20 ± 0.50 




3.30 ± 0.30 


3D Procyon model 


0.95 ± 0.25 


1.45 ±0.25 


3.10 ±0.40 


4.45 ± 0.20 


2.30 ± 0.30 


3.30 ±0.15 


(d3t65g40mm00n01) 















Notes: a:'Holweg er et alj il978l) . full-disk values interpolated: f-(full-disk)=^2(disk-center)/2 + f-(limb)/2; 
b: lsreffen 1 1985.) 



References 

Asplund, M., Nordlund, A., Trampedach, R., 

AUende Prieto, C, & Stein, R. F. 2000, 

A&A, 359, 729 
Barklem, P S., Anstee, S. D., & O'Mara, B. J. 

1998, PASA, 15, 336 
Freytag, B., Steffen, M., & Dorch, B. 2002, 

AN, 323,213 
Holweger, H., Gehlsen, M., & Ruland, F 1978, 

A&A, 70, 537 
Melendez, J., & Barbuy, B. 2009, A&A, 497, 

611 

Steffen, M. 1985, A&AS, 59, 403 
Wedemeyer, S., Freytag, B., Steffen, M., 

Ludwig, H.-G., & Holweger, H. 2004, A&A, 

414, 1121 



